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Abstract 

' In a complex scattering system with few open channels, say a quantum dot with leads, the correlation 

' properties of the poles of the scattering matrix are most directly related to the internal dynamics of the 

system. We may ask how to extract these properties from an analysis of cross sections. In general this is 
^ ' very difficult, if we leave the domain of isolated resonances. We propose to consider the cross correlation 

D I function of two different elastic or total cross sections. For these we can show numerically and to some extent 

also analytically a significant dependence on the correlations between the scattering poles. The difference 
between uncorrelated and strongly correlated poles is clearly visible, even for strongly overlapping resonances. 
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Q> ■ 1 Introduction 

. Starting from Bohr's compound nucleus suggestion, the idea of considering the dynamics in the interaction 
' region almost separately from the scattering process has been proven very successful in different fields. Wigner's 
R-matrix method iQ gives the formal background to the separation of internal dynamics and "free" motion in 
the channel space. Based on this idea we study whether chaoticity or integrability of the internal dynamics can 
be detected in the scattering data. For this purpose we apply Fourier transform techniques, which have proved 
' successful in spectral analysis [|], |[ to total and partial cross sections. 

To build the scattering ensembles, we shall assume that there are no correlations between channel space and 
the internal structure. This assumption is usually well fulfilled for systems with topological chaos, while it is 
often not fulfilled for integrable systems Nevertheless, we use this assumption in order to compare the chaotic 
and regular case in a direct and minimally biased way. Any differences we then find are basis independent and 
^ , minimal in the sense that correlations would usually increase the dissimilarity to the chaotic case. Therefore, 
we use orthogonally invariant random matrix models to describe the internal structure. For the chaotic case the 
choice is obviously the Gaussian orthogonal ensemble (GOE) if time reversal symmetry is conserved. Following 
' Berry and Tabor Q we associate integrability with a random Poissonian spectrum, thus excluding harmonic 
oscillators explicitly. For this case the Poisson orthogonal ensemble (POE) was proposed some years ago 0. 
The case of time reversal symmetry breaking can be treated analogously using the unitary ensembles, but this 
will not be discussed in the present paper. 

We consider three different absorption regimes. For weak coupling the resonances are isolated, and conven- 
tional spectral analysis is satisfactory; if we wish we may add an analysis of the widths. Then follows what is 
usually called strong absorption, where we find overlapping resonances, but where the individual transmissions 
from all channels are considerably smaller than one. Next we have the case where the transmissions are close to 
one. To differentiate the latter two regimes, we shall speak of strong transmission in this case. It corresponds to 
the semiclassical limit, where tunneling effects become negligible. 

For level spectra as well as for intensity spectra, Fourier transform methods and the so called "correlation 
hole" have been very successful tools to identify the effects of integrability and chaos |^ ^, ^, [sj. Yet it is 
not clear to what extent such an analysis remains useful, when the resonances are no longer isolated. An 
exact theory is only available for two-point functions of S-matrix elements in the GOE case. This so called 
Verbaarschot-Weidenmiiller-Zirnbauer (VWZ) integral [|j allows to calculate correlation functions of total cross 
sections by means of the optical theorem. In contrast to that, we have no such theory in the POE case or for 
correlation functions of partial cross sections. In these cases, one had to fall back on the so called "Breit-Wigner" 
approximation px| ], which becomes valid in the weak coupling limit. To extend the validity of this approximation 
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we take advantage of the well known fact that the Satchler transmission matrix ||Tl| or, in absence of direct 
reactions as in our case, the transmission coefficients are the only way in which the coupling strength influences 
the physically relevant quantities. By using the transmission coefficients rather than the coupling constants as 
input, which amounts to a rescaling, we are able to extend the validity of this approximation to the regime of 
strong absorption, where the resonances are overlapping, but the absorption in each channel is still moderate. 

We perform numerical simulations with two purposes: to check the range of validity of the rescaled Breit- 
Wigner formalism, and to find situations, which show significant differences between GOE and POE. We will 
concentrate on Fourier transforms of auto and cross correlation functions, and we shall see that in particular for 
the latter the differences are in some cases very promising. In particular we find that cross correlations between 
cross sections from different channels show strong signatures. 

In Sec. ^ we present the model we use. In Sec. ^ we introduce the correlation functions of S-matrix elements 
and cross sections, our basic tool for the statistical analysis of the scattering systems. Next we discuss the 
"rescaled Breit-Wigner" approximation in Sec. ^. It allows to get results in closed form for the correlation 
functions between total or partial cross sections. In Sec. ^ we derive a twofold integral expression for the Fourier 
transform of the VWZ integral |9). This simplifies the numerical treatment considerably, and it is used in Sec. |[ 
to test the validity of the rescaled Breit-Wigner approximation. After these theoretical considerations, we turn 
to the numerical study of the two scattering ensembles, the POE and the GOE. This is done in Sec. 0, which 
is divided into two subsections. The first deals with correlation functions between total cross sections, and the 
second with correlation functions between partial cross sections. Section H contains a short summary. 



2 The scattering ensembles 

We wish to construct scattering ensembles for the two contrary cases, where the dynamics in the interaction region 
is predominantly integrable, or completely chaotic. We do this under the assumption that the scattering system 
may be separated into an internal part restricted to a finite interaction region and an external part described by 
some superintegrable Hamiltonian. Though some complexity may show up in the coupling of the two parts also, 
i.e. in the mismatch of channel functions and internal functions, the dominant part of the complexity should be 
contained in the internal part. Then the subsystem that describes the complex internal dynamics has a discrete 
spectrum, such that its statistical properties can be modeled with an appropriate random matrix ensemble. In 
order to construct the scattering ensembles, we fix the external part and use standard techniques, originally 
introduced to describe nuclear compound reactions |12, |lj|, to assign to each element from the random 



matrix ensemble a scattering matrix. Thus we obtain a set of scattering matrices provided with the measure 
inherited from the original random matrix ensemble. For simplicity we will denote the scattering ensembles 
obtained from the GOE and the POE by the same names whenever there is no danger for confusion. 

To be more precise, we consider complex scattering systems with many, possibly overlapping, resonances, 
where the S-matrix can be cast into the following form: 

S{E)^l^iV^^^V , H^H,^,-'-VVK (1) 

Here i?i„t is a real, symmetric N x N matrix which describes the internal dynamics, and is a real N x M 
matrix, describing the coupling to the M channels. The matrix H is the so called effective Hamiltonian |l^. In 
order to arrive at Eq. (Q) it is assumed that the coupling matrix elements between channel states and internal 
eigenstates are energy independent |Q. Furthermore one should either neglect the direct reactions or perform 
an Engelbrecht-Weidenmiiller transformation [|l^ if it is necessary to take them into account. The effective 
Hamiltonian H can be diagonalized, such that its eigenvalues Ej = Ej — irj/2 give the positions and widths of 
the resonances if they are isolated. In the eigenbasis of H, the S-matrix elements can be written as 

^ f/ V 

SatiE) = <5.,-i5:^^, (2) 

T E — Ej 

V = A^V , A'^HA = diag (e 



This equation shows that the complex poles of the S-matrix are precisely the N eigenvalues of H. If the coupling 
matrix elements are small enough, their real parts are well approximated by the discrete levels of -ffi„t, whereas 
their imaginary parts are given by the diagonal elements of VV^ . This amounts to the Breit-Wigner approxi- 
mation, which results from applying first order perturbation theory to the effective Hamiltonian H. 
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In this paper we consider two scattering ensembles: the GOE and the POE. Both are invariant with respect 
to orthogonal transformations. Hence in the eigenbasis of H^^^ the M channel vectors are random orthogonal 
vectors. In practice, we use independent random vectors with Gaussian distributed components for V, which are 
orthogonal only up to order 0{N~^). However, as we used relatively large matrices N = 300, the violation of 
the orthogonal invariance had certainly no noticeable effect on our numerical results. Hence, for both ensembles 
the nonzero eigenvalues of VV'' (given by the norm squared of the column vectors), and the level density of iJi„t 
are the only independent parameters. 

In the GOE case, the elements of the diagonal matrix i?i„t are distributed according to the joint probability 
distribution of the GOE spectrum so that for large N, the level density approaches the semicircle distri- 
bution. In the POE case, the elements are independently distributed, and in principle the level density can 
be of any form. Our main objective is the distinction between both ensembles, based on the observation of 
correlations. Hence we find it convenient to use the semicircle distribution in the POE case also. 



The main theoretical tool for the calculation of correlation functions will be the rescaled Breit-Wigner ap- 
proximation, introduced in Sec. ^. In fact, it can be applied in much more general situations. This is of particular 
importance for partially integrable scattering systems, where the assumption of orthogonal invariance does often 
not hold. Then the distribution of the matrix elements of V are typically very different from simple uncorrelated 
Gaussian distributions. 



In cases as they are studied here, the openness of the scattering system is commonly described, borrowing the 
terminology from the so called "optical model" (c/. and references therein), which was originally developed 
to describe nuclear compound reactions with two well separated time scales. Consider the partial cross section 
in appropriate units, which is given by: aab = l^ab — Sab\'^- Then one defines the optical partial cross section as 
^ab — \^ab — ('S'ab)P, where the different time scales are used to obtain a well defined average S-matrix {Sab) by 
averaging over an appropriately chosen energy window. Here (Sab) is simply defined by the ensemble average, 
avoiding in this way any arbitrariness. The openness of the scattering system is then characterized by so called 
"transmission coefhcients," defined for each entrance channel a: 

M 

Ta = {at} ) ~ '^i:! - 1 - E I (^-) I' ' (3) 

where ata} = X^fli ^ac is the total cross section, and a^^} = J2cLi ""ac' i^ ^he total optical cross section, with 
respect to the entrance channel a. The unitarity of the S-matrix leads to: 

ai:} ^2{l^Re Saa) , (4) 

which is sometimes called the "optical theorem." 

The scattering ensembles defined above, i.e. the GOE and the POE, are completely characterized by the 
average level distance d = (Np)~^ in the center of the spectrum of -ffi„t, and the variance of the coupling matrix 
elements (V^), which are independent of i due to orthogonal invariance. From those we define the following 
dimensionless coupling parameters: 

^'^ = iiV'a)- (5) 

Here and in what follows, the brackets (. . .) stand for the ensemble average. As discussed in the Appendix, the 
average S-matrix is diagonal. Its diagonal elements, the transmission coefhcients, and the coupling parameters 
are related to each other: 

which is, however, true in the center of the spectrum only (a more detailed discussion is given in the Appendix). 

As mentioned in the introduction, we distinguish three different regimes. Now, these can be defined more 
quantitatively in terms of the transmission coefficients: The first is the weak absorption regime, where the 
resonances are still well separated, so that -^q ^ 1- N^ext comes the regime of strong absorption, where 

the resonances overlap, but the transmission from each channel is still small: X^fli > f , Va : Ta <C 1. Note, 
that this implies in particular, that the number of channels M is large. Finally we have the regime of strong 
transmission: Va : < I where the transmissions in all channels are close to one. 
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3 Correlation functions of total and partial cross sections 



Correlation functions arc our principal tool for the statistical analysis of total and partial cross sections. We 
distinguish between autocorrelation functions, where one cross section is correlated with itself, and cross correla- 
tion functions, where two different cross sections are correlated with each other. First we define the correlation 
functions in general, in order to introduce our notation. Then we use the optical theorem, to relate the corre- 
lation function of two total cross sections to one of corresponding S-matrix elements. For the latter, the VWZ 
integral provides the exact result in the GOE case. In the POE case an exact result exists only in the one 
channel case and we use the rescaled Breit-Wigner approximation there. For partial cross sections, no exact 
theory exists at all. In this case, we first use the so called "diagonal approximation" to express the partial cross 
sections in such a form, that the rescaled Breit-Wigner approximation can be applied. 

Dealing with matrix ensembles, it is convenient to define the correlation functions as ensemble averages, 
rather then energy averages. Therefore, we eventually have to face the ergodicity question |]l9|| , which is unclear 
in the POE case. Note, however, that in quantum dot experiments ensemble averages may actually be the 
relevant ones We calculate the correlation functions always in the center of the spectrum (see Appendix), 
where we set E — 0. Given then two complex functions / and g of the energy, we define the correlation function 
as follows: 

C[/,5]M = {.f{^) g(f )) - (/(-f^)) )) . (7) 

Here d is the mean level distance in the center of the spectrum of i?i„t which is assumed to be constant on the scale 
where we expect correlations. Note that, there is no unfolding involved. The mean level distance d simply serves 
as a convenient energy scale. For the discrete spectrum of some random Hamiltonian: f{E) = ~ the 

autocorrelation function becomes: C[f, = 1 + S{llj) — i^2(w), where I2 is the two-point cluster function as 



defined in Ref. 1 17 



We will mainly analyze the correlations in the time domain, and by consequence deal with the Fourier 
transforms of correlation functions. In general, we denote the Fourier transform of a given function of the energy 
f{E) by: 



/(t)-.F[/](t)=y d^e^--'/M), (8) 

where the Fourier transform is taken with respect to the energy measured in units of d, and the factor 27r in the 
exponent assures proper normalization. For two spectral functions / and g, the following relation holds: 

C[f,gm = i {{f{-t) g'it)) - ifi-t)) (g'it))] , (9) 

where the functions /'(x) and g'{x) are equal to f{x) and g{x) inside the interval \x/d\ < L/2 and zero outside, 
and the limits N, L ^ 00 , L/N — > are taken. Equation is based on the convolution theorem 21| applied to 
the fluctuating parts of the spectral functions / and g, where the convolution integral is expressed as a correlation 
function as in Eq. (|^) assuming stationarity. The limit A'^, L — > 00 is necessary, to allow the correlation function 
to go to zero quickly enough, so that the Fourier integral of the correlation function is well defined. The 
limit L/N serves to obtain stationarity in the interval where the correlation function is calculated. In 
particular, the average level (or resonance) density and the average S-matrix should not vary noticeably in this 
interval. Equation (^ is used in the numerical calculation of the correlation functions. It turned out, that 
L = N/2 , N = 300 already gives well converged results. 

Note that we measure the energy in units of the mean level spacing d. As a result, the argument of a 
correlation function is dimensionless, and so is the argument of its Fourier transform. Nevertheless the latter is 
denoted by i, as its significance is still time — though measured in units of d^^. 



Total cross sections Consider the correlation function of two total cross sections fitot and atol with possibly 
different entrance channels a and b. As the total cross sections depend linearly on the respective diagonal 
S-matrix elements [see Eq. (^] the correlation function can be expressed as follows: 

C^Wtat^^tatl = C[Saa + S*a, Sbb + S^b] 

= 2Re {C[Saa, Sbb] + C[Saa, Sbbl) 

= 2ReC[Saa,S;b]- (10) 
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The correlation function of nonconjugated elements C[Saa, S^b] vanishes |2^. Relation ( pO| ) is essential, as it 
relates experimentally accessible quantities to analytical results For the Fourier transform of Eq. (10) 

we get: 

(?K.,^L]W = 2J^ReC[Saa,St,]{t) 

= C[Saa, S^b]{t) + C[Saa, 'S'bf,]( — 

= C[Saa,S*Jit) , (11) 

where it is assumed that t > 0. Then C[Saa, Slij]{~t) vanishes, because of its negative argument p4] (see also 



Sec. 4.1). For the sake of brevity let us think of a correlation function and its Fourier transform as a single object 
represented in the energy domain or in the time domain, respectively, and call it simply "correlation function" 
or "C-function" in either case. 

Partial cross sections Partial cross sections are given by aab = \Sab — Sab\'^- The theoretical treatment 
of correlation functions of partial cross sections is complicated by the fact that one has to average over a 
product of four S-matrix elements. The insertion of the S-matrix elements as given in Eq. (^ leads to a double 
sum of resonance terms. As an exact analytical treatment seems to be impossible, we employ the diagonal 
approximation, which consists in retaining the diagonal terms of the double sum only. This is justified for 
sufficiently weak coupling and leads to: 

N ^2 

The r.h.s. can be written as the imaginary part of a function jab{E), which has the same pole structure as the 
S-matrix: 

cj',,{E) = -2Im/,fc(i?), 

Due to the linear relation between the diagonal approximation cr^^ and the spectral function jab, we may again 
express correlation functions of the former by corresponding correlation functions of the latter. In fact, we have: 

C[aab, cj,d]{Lo) « CKfc, - 2KeC[!ab, Ui^) ■ (14) 

It remains to calculate the correlation function of Jab and J*^. This will be done in the following section using 
the rescaled Breit-Wigner approximation. In order to calculate the Fourier transform of Eq. (p^), we note that 
the pole structure of Jab{E) is the same as that of the true S-matrix, so that ^[/at, /*J(— t) is again zero (for 
< > 0). Therefore, we obtain: 

C[aab, <7cd] (t) « CK„ a'J (t) CUab, U it) • (15) 



4 Rescaled Breit-Wigner approximation 

In this section we calculate the correlation function of two arbitrary S-matrix elements using the Breit-Wigner 
approximation, followed by a phenomenological rescaling procedure. To this end Eq. (^ is used, which means 
that we first calculate the Fourier transform of the respective S-matrix elements (this can be done exactly) , and 
then we average over the resonance parameters {Ej, Tj, Vji, . . . , VjM}i<j<N ■ The average over the real parts of 
the S-matrix poles {Ej} can still be done in a formally exact manner, but then we have to use the approximation 
mentioned above, in order to finish the task. To obtain the correlation functions of partial cross sections, the 
same steps have to be done with the matrix elements jab instead [see Eq. ([T^)]. 
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4.1 Formally exact treatment 

Here, we do all those steps of the calculation which are exact. We first calculate the Fourier transform of one 
S- matrix element, using its pole expansion (^). For i > 0, we get: 

— 1 _ _ g-27ria;t 



-i A /"^/^ 

Sab{-t) = -ry.^i-v,b / dco 

" ~1 J-L/2 



= ^X^VSaVS-6e-2-^^*/rfe-r.*/^ (16) 

Here it was used, that the poles with real parts outside the integration region do not contribute to the Fourier 
transform, and for those inside, it is well justified to extend the integration up to infinity because Fj/L <C f . 

In the same way, we may obtain an analogous expression for Scd{t)- The Fourier transforms of the average 
S-matrix elements are taken into account in the numerical calculation only. Here by contrast, we notice that 
the average S-matrix elements are almost constant in the integration interval, which means that in the limit 
TV, L — > cxD their Fourier transforms will become ^-functions situated ai t = 0. Thus they play no role in the 
current calculation which is restricted to t > 0. Note that we may define C[Sab, S*^]{0) = limf^o C[Sab, S*^]{t). 
Inserting the expressions for Sab{—t) and S*j^{t) into Eq. and ignoring the Fourier transforms of the average 
S-matrix elements, we obtain: 



1 ^ 

C[Sab,S*J{t) = J2v,aV,tV:,V:,e 



~2m(Ej~Ek)t/d „-Tr(rj+rfc)t/rf 

= {{y^aVo^y;cV;d e-'^^^^/^) + {L-1) {v,aV,,V,*A*d e-^(^^+^^)*/'' ,^2.nE,^EMd^^ . 

(17) 

As the ensemble average is invariant for any permutation of the resonance indices, the double sum can be 
evaluated. In the final expression the resonance indices j ^ k are arbitrary. Note that C[Sab, S*j\{t) vanishes for 
t < 0, because in this case both Fourier transforms Sab{—t) and 5*^(t) vanish, as can be easily seen by applying 
the residue theorem (cf. also p^). 

At last, we average formally over the real parts {Ej} of the S-matrix poles. For fixed values of the partial 
amplitudes and the total widths, the average over —L exp[— 27ri(ii'j — Ek)t/d] for L 00 gives the two-point 
form factor ||l^ of the random sequence {Ej}, which we denote by 62 (i)- In general, 62 (i) still depends on the 
parameters fixed. In the weak coupling limit, however, the positions of the resonances on the one hand, and the 
partial amplitudes and total widths on the other hand, become statistically independent, so that 62(0 converges 
to the two-point form factor of the closed system. After all we may write: 

C[Sab, S*J{t) = 1 {{v,aV,,V;,V;, e-^-r.*/^^ _ {v,aV,bVC,Kd e--^^^+^^^'/'' Ut))] . (18) 

This is so far an exact but rather formal result. However, it clearly shows that the correlation function is no 
direct measure for spectral correlations. The first term in Eq. (|l8|), which may very well dominate the correlation 
function, contains the parameters of only one single resonance. Therefore, it cannot describe correlations between 
different resonances that are the ones we are really interested in. It is typically a monotonously decreasing 
function, where the decay is governed by the average width of the resonances. It is the second term in Eq. (p^, 
which contains the parameters of two different resonances. It vanishes completely if there are no correlations 
between them. For our investigations it is important to find situations, where the first term is relatively small, 
so that one may retrieve as much information as possible on the correlations between different resonances. 

4.2 Approximation 

In order to evaluate the remaining averages in Eq. we have to introduce some approximations. To this end, 
consider the weak coupling limit: Va : Ka — > 0. Using first order perturbation theory in the expression for the 
S-matrix, Eqs. (|l|) and (|^), we get the pure Breit-Wigner approximation for the S-matrix: 
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This amounts to make the foUowing replacements in the pole expansion of the S-matrix (^: 

M 

Vja^V.a, r,-^V^?, E,^e,, (20) 

where Sj are the eigenvalues of -ffi„f Hence, in order to obtain the correlation functions in the Breit-Wigner 
approximation, we simply do the same replacements in Eq. (|l^). Then the partial amplitudes become real 
uncorrelated Gaussian random variables, the total widths become simple functions of the partial amplitudes, 
and 62(0 becomes the two-point form factor 62 (^) of the spectrum of i/i„t: 

C[Sa,,S*Jit) « ^ {{v,aV,,V,,V,, e-^-^^K ^/'^) _ {v,aV,,VucVu, e--^.(K+yL)^/<^)b^{t)] . (21) 

The remaining Gaussian averages are relatively simple, so that in many cases the respective correlation function 
can be calculated in closed form. Note that the averages are different from zero, only if all partial amplitudes 
appear in even powers. 

Unfortunately, the pure Breit-Wigner approximation drifts quickly away from the exact result, as the coupling 
to the continuum increases. The following phenomenological procedure improves the approximation considerably. 
It consists in rescaling the variance of the partial amplitudes as follows: 

TtV^ • (22) 

We call the result the "rescaled Breit-Wigner" approximation. As shown in Ref . |^ , it leads to partial fluctuating 
cross sections of the Hauser-Feshbach type ||2q|, showing elastic enhancements of 2 (GOE case) and 3 (POE case) 
in agreement with earlier theoretical resultepst |2^ , |29|| . The occurrence of Ericson fluctuations |Q with the 
correlation length Tq = ^X^^i '7o/(27r) is also correctly described. In Sec. Ijit is shown that in the time domain 
this approximation gives results for the correlation functions between S-matrix elements which become exact as 
t 0. 

Note that the rescaled Breit-Wigner approximation can be applied in a wide range of situations, while 
it conserves the simplicity of the pure Breit-Wigner approximation. This makes it a valuable tool for the 
statistical description of complex scattering systems. This approximation can be justified to some extent with 
the following reasoning: It is well known, that the properties of a scattering system are determined by the Satchler 
transmission matrix [ pT| , if the entire process is occurring on two different time scales, one associated with direct 
processes and one with long time processes also called compound processes in Nuclear physics. In the absence 
of direct reactions or after an Engelbrecht-Weidenmiiller transformation , this implies the dependence on the 
transmission coefhcients only. They are directly related to the variances of the partial amplitudes as mentioned 
above. The rescaled approximation thus implies that we use the transmission coefficients of the system rather 
than the coupling constants. This can be viewed as a nonperturbative input, or as using phenomenological 
parameters. 

Total cross sections 



According to Eq. (|1 1[) , the C-function between two total cross sections is equal to the C-function between the 
respective diagonal S-matrix elements. Hence C[a['^} , ](i) is given by the r.h.s. of Eq. (18) setting b = a and 



d = c = b. In order to calculate this C-function in the rescaled Breit-Wigner approximation, we use the Eqs. d) 
and (|^) to express the ensemble averages as integrals over normalized squared amplitudes gi, ■ ■ ■ ,gM'- 

C[al:l<^itm « {{9a9t e-^ *) - (.9. e^^ ^Z^) {g, e'^ */^) b,it)} , (23) 

with G = X^fli 3c- The normalized squared amplitudes gc are distributed as the random variables Vj^/ {V^), 
i.e. they are Porter-Thomas distributed |^ : 

p(g^) ^ ^= e-3^^\ c=l,...,M. (24) 

Now, the remaining averages in Eq. (|23| ) can be calculated easily |2^ . For the auto C-function of the total cross 
section we get: 

C[ath{t) ^T^3i^ + 2r„t)-5/2 + 2r,t)-V2 _ (1 + T^ty^ + T,t)-' b,it) I . (25) 
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Note, that this result is very similar to the result of Ref. |Q for the spectral autocorrelation function of the 
photodissociation cross section in weakly coupled chaotic systems (there the rescaling is used also, though 
without mentioning it; note that these authors obtain a rescaled formula for their specific case directly as an 
approximation to exact results ||3^). For the cross C-function of two different total cross sections, we get: 

C[at\^^'}m ~ TaT, \ (1 + 2T„t)-3/2(l + 2^1)-^'^ [] + '^T,t)-^/^ 

-(1 + T„i)-2(i + r,i)-2 [] (1 + T,ty' 62 W \ ■ (26) 

In what follows, we assume that all coupling strengths are equal: Va : k = Ka, T = Ta- For convenience, the 
numerical analysis in Sec. |^ is restricted to this case only. The Eqs. ( p5| ) and (26) simplify considerably and can 
be combined into a single one: 

C[at\'ji'J] « T2 {(1 + 2Sa,) (1 + 2Tt)-'-''/' - (1 + Tty'-^' b,{t)} . (27) 

To lowest order in the transmission coefficient T the r.h.s. becomes [1 + 25 ab — &2(^)] m agreement with 
results on intensity weighted stick spectra ^ and the asymptotic behavior of the exact analytical result for 
the GOE case (see Eq. ( |59| ) in Sec. ||). The difference in the C-function at small times between GOE and FOE 
is known as the correlation hole. In order to quantify it, we use its size at i = relative to the maximal size of 
the C-function in the FOE case. From Eq. ( ^ ) it follows that the correlation hole is 1/3 in the case of the auto 
C-function (a = 6), while it is one in the case of the cross C-function (a ^ b). 



Partial cross sections 

According to Eq. (|l5|) , the C-function of two partial cross sections in the diagonal approximation is equal to the 
C-function of the matrix elements Jab and /*^, defined in Eq. (p^). Following the same lines as in the case of the 
true S-matrix, we arrive at the following expression for the latter C-function: 

C[U,a{t) « ^ e--^^ ^Z'^) - e- (^^+^^) ^^(O)} , (28) 

which is the analog of Eq. (|l^). To this expression we can apply the rescaled Breit-Wigner approximation. 
The partial widths jja, are replaced by Porter-Thomas distributed random variables with average value (7^0) = 
Ta (i/(27r), the total widths become the sums of the partial widths, and the two-point form factor 62(0 becomes 
the two-point form factor of the spectrum of i/i„t. Thus we get: 

CKfc,a,,](0«T,nTO{(^^i^|^e"«*)--(^e-«*/2)(^e-«*/2) b^it)} , (29) 
where G ~ X^c^c^c- The ensemble average can be performed, using the following identities: 

/ da' e-"'« , 

/ da' / da" c-°'"^ , a > . (30) 

J a J a' 

Exchanging the integration on a' and a" with the ensemble average, we can do the ensemble average analytically, 
and we are only left with the integrals over the auxiliary variables a' and a". With the exception of some special 
cases, those integrals have to be evaluated numerically (for details, see Ref. [^). In the case of equal coupling 
strengths, however, we obtain the following analytical expression: 

<?K>> ^.^1 w - { (1 , i).;,T/4/M) ' " + - ° " + ■ 

Here, A and B depend on the actual combinations of partial widths in Eq. (|29|), i.e. whether any of the channels 
involved coincide or not. The following cases occur in the present paper: 
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• The autocorrelation function: 

CKfc] : A = 9 + 965ab, B ^l + 8Sab. 

• The cross correlation function between two different elastic cross sections: 
C[aaa,<Jbb] : A = 9, 5 = 9. 

• The cross correlation function between an elastic and an inelastic cross section: 
C[<Jaa,aab] : A - 15 , S = 3 . 



According to Eq. (|31|), the coefficients A and B determine the relative depth of the correlation hole in the limit 
of many channels M -> oo. The prefactor (1 + 2/M) + 6/M)(l + 4/M)] in front of A may, however, lead to 
a considerably deeper correlation hole, as long as M is not too large. 

For M ~> oo and MT fixed, the algebraic decay in Eq. (|3l]) turns into an exponential one: 



(1 + 2Tt)"2"*^/2 _ (l+Tt)-^-'^^ - exp(-A/rt) . (32) 
This means, that we obtain Ericson fluctuations l30| in this limit. The reason for the occurrence of Ericson 



fluctuations can be understood from Eq. (29): Only in the case of many channels and small transmission 
coefficients, does the central limit theorem lead to negligible fluctuations of the total width around its average 
value. Then we may treat C? as a constant, which leads immediately to the expected exponential decay. Note 
that this implies that Ericson fluctuations are no reliable signature for chaotic scattering, because the central 
limit theorem may work, even if the partial amplitudes are not Gaussian distributed, thus leading again to 
Ericson fluctuations. In such a situation, the deviation from the exponential decay due to b2{t), would be the 
only reliable signature of the chaotic dynamics. 



5 The Fourier transform of the VWZ integral 

In this section we derive a general formula for the Fourier transform of the VWZ integral 1^. This allows to 
obtain exact results for the correlation functions of total cross sections in the GOE case. We need these results, 
in order to check the accuracy of the rescaled Breit-Wigner approximation. Apart from that it will turn out, 
that numerically it is much easier to calculate the Fourier transform than the original VWZ integral: 



8 J Jo Jo ^Xiil + Ai) A2(l + A2) (A + Ai)2(A + A2) 



2 



X 



Ilab = TaTb 



Y[ I, rr ^ rr T=\ {^ab^cd Aa Ac + {^ac^bd + ^ad^b,^ Ilab} , (33) 

e=i V (1 + JeAl)(l + TeA2) 

Ai A2 2A 



l + r,Ai l + TaA2 1-T„A^ 
Ai(l + Ai) , A2(l + A2) , 2A(1-A) 



(1 + TaAi) (l + TbAi) (l + TaA2) (l+TbA2) (1 - T,A) (1 - nA) 

The main interest in the VWZ integral has been the calculation of average fluctuating cross sections ^ 
p6{ , which corresponds to the case w = in Eq. (|3^) . Only a few papers treat the w-dependence of the VWZ 
integral [ p4| |3^ , |38[ | , and even then the analysis was usually restricted to particular limits such as w — > (t ^ cxj 
in the time domain), or many channels: M ^ oo and Ericson fluctuations. Quite often the existence of the 
correlation hole was simply ignored. 

In the present analysis, the correlation hole and its dependence on the coupling strengths is of great impor- 
tance. It is needed to distinct regular from chaotic dynamics. Therefore, we will analyze the Fourier transform 
of th e VW Z integral in some det ail. We will also observe the behavior of the correlation hole in the limits: i ^ 



(Sec. |5j) and Va : Ta ^ (Sec. |5J). 



Let us start our derivation by applying the Fourier transform (g) to Eq. (p3|). We exchange the Fourier 
integration on lo with the integrals on A, Ai, and A2. For the Fourier integration we then simply need to 
calculate: 

dw exp{27ritj [i- (Ai + A2 + 2A)/2]} 

= (5[t-(Ai+A2 + 2A)/2] . (34) 
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The (5-function can be used to remove the A-integral. To the remaining double integral, the following transfor- 
mations are applied: 

\ \ Ai + A2 

Ai , A2 r = , s = A2 - Ai 

followed by s x = - . (35) 



This leads to: 



where 



1 r* *^ 

C[Sab, S*Jit) = - dr{t- r)(r + 1 - t) J] [1 - T,(t - r)] U{r) , (36) 

4 Jmax(0,t-1) „_i 



U{r) = 2 dx = , (37) 

(t^ - r2 + x)2 ^x{x + 2r+l) ^JUtLli^ + ^^eT + T^x) 

Aa = 2Ta^fr^a( ^±1^ + l—L -^,and (38) 

^ ''\l + Ta{2r + Tax) l-Ta{t-r)J ' ^ ' 



Uab — 2TaTb 



TaTbX^ + [TaTbr + {Ta + Tb){r + 1) - l]x + r{2r + 1) , {t - r)(r + 1 - t) 



(1 + 2Tar + T^x){l + 2nr + T^x) [1 - T,(t - r)] [1 - n{t - r)] 



(39) 

In order to remove the l/\/x-singularity in the integrand of C/(r), we substitute consecutively: 

x = {y-l)b/2 , y = {z + z^^)/2 , z^2u+\ where & = 1 + 2r . (40) 



This gives: 



jj, . . r dw SabScd AaAc + (SacSbd + SadSbc) ^ab hu^ 

U{r) = 4 / -^^^^^^^^^^^^= , X = . (41) 



The Eqs. ( |36| ) and ( [11[ ) together with the Eqs. (38) and ( |39D form our final result for the correlation function 
C[Sab, S*j}{t) (it is understood, to replace x by bu^/{2u + 1) wherever it occurs). One may readily use these 
formulas for numerical calculations. All problematic singularities have been removed from the integration region. 
Our result if followed by a fast Fourier transformation back to the energy domain, may even be a quite efficient 
way of computing the original VWZ integral. 

The Figs. |l| and ^ show examples of the Fourier-transformed auto and cross correlation functions of the total 
cross section, calculated with the help of Eq. (|3^), etc. The coupling strength for all channels was fixed to the 
value K = 0.1, while the number of channels M was varied. The qualitative features of the C-functions are almost 
the same for any M. In the first place, M determines the falloff at large times. Note, however, that in the case 
of extremely small total transmission, we would obtain a "true" correlation hole at t = (i.e. a positive slope). 
See also the related discussion below Eq. (p7|)] which is concerned with the rescaled Breit-Wigner result. In all, 
it permits us to restrict the numerical studies in Sec. ^ to the case M = 2. 

Figure [l] shows the autocorrelation function divided by T^, which starts at t = with the value 2, independent 
of the number of channels. It decays with time, the larger the number of channels, the faster. This is very 
reasonable, as the decay must be related to the total transmission J^a -^f ~ MT. Note that there is no qualitative 
signature (positive slope at i = 0) of the correlation hole. 

In Fig. H the cross correlation function is shown, which starts at t = with zero. It increases linearly with 
time, and it seems that now the slope is independent of the number of channels. At large times, the cross 
C-function must go to zero in the same way as the auto C-function shown in Fig. ^ Hence, after reaching a 
maximum that decreases and moves to the left with increasing M, the cross C-function decays to zero. Therefore, 
we see a clear signature of the correlation hole. 



To conclude this section we will consider two different asymptotic limits: First we will prove, that for t 
the asymptotic behavior of the various C-functions is independent of M. Second we will show with very high 
numerical precision, that for small transmission coefficients, Va : Ta ^ 0, the cross C-functions coincide (if 
properly normalized) with the two-level form factor b2(t) for the GOE. 
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t 

Figure 1: The exact result ( |3^ ) for the autocorrelation function C'[crtot](^) divided by . While k = 0.1 is kept 
fixed, the number of channels is varied: M = 2 (solid line), M = i (long dashed line), M = 10 (short dashed 
line), M = 20 (dotted hue), and M = 40 (long dotted line). 




t 

Figure 2: As in Fig. |^, but for the cross correlation function C[a'^^^, '^tut\{^) divided by T^. 
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5.1 The limit of small times 

In the limit t ^ 0, the Eq. ( ^6|) can be approximated by: 



C[Sab, S'cd] (t) 



J I dg{l-g) U(tg) 



dg (1 - g) U{g) 



(42) 



where U{g) = t^U{tp). Hence setting r = tp m Eq. ( ^ ) the following limits can by taken: h ^ 1, 2w+l'-^ 1, 
and cc ~ w^. Then substituting u = tv we get: 



U{g) - - / dw [SabScd ^a^c 



{5achd + SadSbc) ^ab] 



dabScd AaAc + {6ac6bd + ^ad^bc) ^ab 



To lowest order the Eqs. (^8|) and 



dw 



simplify also: 



(1-^2 



,,2^,2 



2Tav/l-Ta< (1 + Tatv"^) ^ 2Ta^l - Ta t 
2TaTb t [1 + {Ta +Tb- l)tv^] - 2TaTb t . 



(43) 



(44) 
(45) 



This shows, that if the C-function contains the A-term only, we get a linear increase starting at i = 0, whereas 
if it contains the Il-term also, we get a finite value at t = 0, but no valid value for the slope. The slope in the 
first case and the finite value in the second are both determined by the same integral Ii : 



C[Sab,S*j] (t) 



4 SabScd TaTc \/(l - Ta){l - Tc) t + 2 {SacSbd + SadSbc) TaTb 



h , 



dv 



— £|2 + t|2) 



2\2 



f dg{l^g) 
Jo 



2p3 



■ arctan | — 
P 



Q 
2p2 



(46) 
(47) 



where = 1 — (? . Splitting 1\ into two parts and substituting x = arctan(gi/p) in the first one, we can evaluate 
Ii analytically: 



h=- 



7r/2 



dx 



X 



1 + sin X 



g In 2 l-ln2 1 

do — ^ = \ = - 

1 + g 2 2 2 



(48) 



In order to contain the A-term only, the C-function must be of the form: C[SaaT S^^^ with a ^ h. Then its 
asymptotic behavior is: 



C[Saa. Sl^] ^ 2 TaTb ^ {1 - Ta){l - Tb) t . (49) 

If the C-function is of the form C\Sab: S'*^] = C[Sab, Sla\: then it contains the Il-term also. In this case we get: 

C[Sab, Sl^m ^{1 + Sab) TaTb . (50) 

These results confirm our observations in Fig. ^ and Fig. ^j. The cross correlation function of two total cross 
sections vanishes at i = 0. This holds for all values of the transmission coefficient. Moreover, the slope only 
depends on the transmission coefficients of the entrance and exit channels, and thus does not change if the 
number of channels is changed. In contrast to that, for the auto correlation function of a total cross section we 
get: C[a['^}]{Q) — 2Ta- Again this value does not depend on the number of channels. 

Comparing the asymptotic results (^9|) and (^0|) to the rescaled Breit-Wigner approximation (|2^) we find 
agreement in lowest order — for the slope in the case of cross correlation functions, and for the value at i = in 
the case of autocorrelation functions. 
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5.2 The limit of small transmission coefficients 

In the case Va : Ta 0, the Eqs. @ and (|l|) simplify to: 

1 



4 

where 



C[Sab. S:M) ^ 7 / * dr{t-r){r + l-t) U{r) , (51) 

inax(0,i — 1) 



^"^'^ ' 2U+1 (t2 _ ^2 + _52i!_)2 • ^ > 







In the same Hmit, the expression for A^ and the one for 11^^ ( |39| ) become 

(t - r){r + 1 - t) + r{2r + 1) 



Ilafa ~ 2TaTb 

Defining the two integrals: 



Aa - 2r,t (53) 



2-^+ 1 



(54) 



r , 2u + 1 

h ^ du ^ 55 

Jo [(<2-r2)(2u + l) + feu2]2 

f u2 

Jo [(t2-r2)(2w + l) + &w2]2 ' 

the asymptotic limit of ?7(7') can be written as 

U{r) ^ 8 {26abScd TaT, h + {Sachd + SadSbc) TaTb [(r2 +2tr + t- <2)/2 - bh] } . (57) 

In order to keep the discussion short, we shall consider cross C-functions of the type C[Saa, S^i,], a ^ b only. 
Then we only need to calculate l2- Defining =t^ — and = b — p^, I2 can be written as follows: 

1. /"■ du 

h = -—Or, 



2p ^ Jo bu"^ + 2p'^u + p2 

11/ 6r+p2 p 

Or, arctan — arctan — 

2p ^ pR\ pR R 

R'^-p^f br+p'^ p\ r(r + l) 

arctan arctan — ' 



2p3i?3 V pR R) 2t2p2^2 

arctan — + ) ' . (58) 



2p3i?3 (r + l)p 2t2p2^2 

Inserting Eq. (|^ into Eq. (|5l|), we get the following result for C[Saai S^f] 

1 



TaTb 



C[Saa, S;,]it) ~ At' / dr{t-r){r + l-t) I2 , (59) 

v'max(0,t-l) 

where I2 is given in Eq. (|5^). We have evaluated Eq. ( ^9|) numerically and compared to the expected result 
1 — 62(^)5 for the discrete GOE-spectrum. The difference is of the order of the machine precision (sa 10^^^) for 
any value of t. 

6 Test of the rescaled Breit-Wigner approximation 

The aim of the present paper is actually twofold. In the first place, we analyze correlation functions of total and 
partial cross sections, to search for signatures of the correlation hole, in particular in the regime of overlapping 
resonances. In the second place we wish to establish the rescaled Breit-Wigner approximation as a simple and 
general tool for the analysis of correlations in different scattering situations. This section is devoted to the latter, 
whereas the the correlation hole will be studied in Sec. |^. 

We use the rescaled Breit-Wigner approximation in those cases, where no exact theory is available, i.e. for 
C-functions of partial cross sections (where we need the additional diagonal approximation) and for C-functions 
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Figure 3: Difference between the rescaled Breit-Wigner approximation ( p7| ) and the exact resuh (|3^) for the 
autocorrelation function C[(Jtot](0 divided by T^. The couphng strength is kept fixed k = 0.1, and the number 
of channels is varied: M = 3 (solid line), M = 10 (long dashed line), M — 20 (short dashed line), and M = 40 
(dotted hue). 

0.05 I I I I I I I I I 




0.5 1 1.5 2 2.5 3 3.5 4 
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Figure 4: As in Fig. ^, but here, the number of channels M = 2 is kept fixed, and the coupling strength is varied: 
K = 0.01 (sohd hue), k = 0.03 (long dashed line), k = 0.1 (medium dashed line), k = 0.3 (short dashed line), 
and K = 1 (dotted line). 



14 





15 



of total cross sections, in the POE case. In this section we use the C-functions of total cross sections in the GOE 
case (for which we have exact analytic results) to check the validity of the rescaled Breit-Wigner approximation. 



The Figs. ^|-|6|show the difference between the rescaled Breit-Wigner approximation ( |27| ) and the exact result 
36| ) for C-functions of the total cross section in the GOE case. This is done for different numbers of channels 
(Fig. H and Fig. and for different coupling strengths (Fig. ^ and Fig. |^). As the error of the rescaled Breit- 



Wigner approximation is zero for t = (see Sec. 5.1), all curves for the auto and cross C-functions begin at zero. 
In the case of the auto C-function the rescaled Breit-Wigner approximation mostly underestimates the exact 
result. In particular at small times, we get maximal deviations, which decrease again as t — > oo. In the case of 
cross C-functions the error curves depend in a more irregular way on the number of channels or on the coupling 
strength. Though the error is of comparable size in all cases, the error curves at t — > seem to behave linearly 
in the case of the auto C-functions and quadratically in the case of the cross C-functions (in agreement with 



our theoretical expectations; see Sec. 5T). Yet the cross C-function vanishes for k — I, while the approximation 
does not and must, therefore, be considered with great reservations. In general, the error is rather an absolute 
then a relative one. We may say that the rescaled Breit-Wigner approximation works quite well up to k = 0.1, 
and that for larger k the approximation still behaves qualitatively well except for the cross C-function for k — 1. 

In Fig. 0, we compare the rescaled and the ordinary Breit-Wigner approximation. The latter is obtained by 
undoing the rescaling and replacing T by 4k in Eq. (|2^). As we have no exact theory in the POE case [Fig. 
0(b)], we use numerical Monte Carlo data instead (c/. Sec. 0). The data may be considered as exact up to 
the small fluctuations, seen in the data points. Though the resonances are well separated for k = 0.1 which 
is used here, the pure Breit-Wigner approximation is already far off the exact result, in particular at small 
times. This is true for the GOE case shown in Fig. 0(a), as well as for the POE case in Fig. 0(b). In fact, the 
rescaled Breit-Wigner prediction in the GOE case is much closer to the POE case, which makes it practically 
impossible to distinguish between correlated and uncorrelated spectra. Note that the auto C-function shown 
here, is monotonously decreasing without showing any qualitative sign of the correlation hole. 

The rescaled Breit-Wigner approximation by contrast, reproduces the exact values of the C-function at t — 
in the GOE and apparently also in the POE case. This allows to detect the correlation hole if present, at least in 
principle (absolute cross sections must be available with sufhcient accuracy). Though the rescaled Breit-Wigner 
approximation underestimates the auto C-function at intermediate times in the GOE case [Fig. 0(a)], it agrees 
very nicely with the numerical data in the POE case [Fig. 0(b)]. 

It seems that the rescaled Breit-Wigner accounts very well for uncorrelated resonances that diffuse inde- 
pendently into the complex plane (similar to what one would expect from an effective mean field theory). For 
resonances which are strongly correlated from the very beginning, the dynamics of the resonances is more com- 
plicated, and the rescaling procedure cannot fully account for that. However, in cases where deviations of a few 
percents (of the maximal C-function value) are still acceptable, the rescaled Breit-Wigner approximation can 
be applied up to quite large coupling strengths, which reach well into the strong absorption regime. While it 
is convenient to consider the absolute error at small times (in the correlation hole region), the relative error is 
more appropriate at large times, i.e. in the tails of the C-function. The latter apparently scales with the total 
transmission X]a=i -^a have checked this separately), in line with earlier theoretical work [Q. 



7 Numerical analysis: GOE versus POE 

Here we perform Monte Carlo simulations for correlation functions of partial and total cross sections. The aim 
is to distinguish the situations, where the internal Hamiltonian has a purely random spectrum (POE), from 
those where it has a spectrum with correlations (GOE). We restrict the analysis to coupling strengths < k < 1 
that covers the whole range for the transmission coefficient. As discussed among St others in Refs. § ||, |o) the 
scattering ensemble may show quite a different behavior when k ^ 1. The results are presented for two open 
channels. Calculations for larger channel numbers have been performed, but they show no significantly different 
behavior. In the case of total cross sections and GOE, this can be seen from the results for the VWZ integral 
obtained in Sec. V (see for example Fig. 0). In the other cases, i.e. POE or partial cross sections, the rescaled 
Breit-Wigner approximation gives the same answer [see Eq. (^)]. In all cases the number of channels merely 
enters as a scale factor. 

The numerical data are obtained using Eq. (0), i.e. we calculate the Fourier transforms of the two cross 
sections in question and then take the ensemble average (400 runs) over their (complex conjugated) product. 
Finally the product of the average cross sections is subtracted to eliminate the peak at i = 0. The dimension of 
the effective Hamiltonian H [see Eq. (0)] is = 300. Where the resulting data are smoothed, this is done using 
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Figure 7: The autocorrelation function C[atat]{t) divided by for M = 2 and k = 0.1. The long dashed 
line is the rescaled Breit-Wigner approximation (|2^), and the short dashed line is the ordinary Breit-Wigner 
approximation (details see text). In (a) we show the GOE and in (b) the POE case. The solid line in (a) shows 
the exact result (|36|), whereas the crosses in (b) show numerical Monte Carlo calculations. 
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Figure 8: The autocorrelation function of the total cross section divided by T^. The GOE versus the POE case, 
for coupling strengths: k = 0.03 (a), k = 0.1 (b), k = 0.3 (c), and k — 1. (d). The circles show the smoothed 
numerical GOE data; the crosses the equally treated POE data. In all cases (a), ... ,(d) the exact VWZ integral 
is used for the GOE theory (solid line) and the rescaled Breit-Wigner approximation for the POE (dashed line). 



a x^-fit to a natural cubic spline. The length of its curve segments was 28 points that corresponds to At « 0.26 
[ pl| . From the spline, 150 equidistant points in t arc finally plotted. Note, that there is no free fit parameter, 
and the curves are not normalized. The factor by which we divide, is directly determined from the input 
parameter k, using Eq. (ra). 



7.1 Total cross sections 

We start the numerical analysis with auto and cross correlation functions of total cross sections. While this case 
is relatively easy to handle in theory, it is usually very difficult to obtain total cross sections from an experiment, 
in particular if two independent total cross sections are needed. 



Autocorrelation functions 

As a starting point consider the correlation hole in the Fourier transform of a discrete energy spectrum 
C{t) ~ 1 — 62(^)1 that distinguishes the GOE from a purely random spectrum (POE). In the limit of vanishing 
coupling K — ^ 0, the auto C-function C[(Jt^t]{t) divided by converges not to C{t), but to C{t) + 2. The reason 
is that the resonances are weighted by Porter-Thomas distributed intensities (or partial widths) |^, which 
leads to a reduced correlation hole of 1/3. 

Figure || shows C[atot\it) /T'^ , for the GOE and the POE case. We find the behavior described above, but 
with the difference, that due to the finite coupling C (t) must be multiplied with roughly an exponential decaying 
function. In the four panels (a)-(d), the coupling parameter is increased from k — 0.03 to k = 1. In the GOE case 
the numerical results are compared with the exact theory, Eq. (^6[), but in the POE case we are dependent on the 
rescaled Breit-Wigner approximation, Eq. (p7|). We find perfect agreement for the GOE which shows that the 
numerical Monte-Carlo calculation is reliable. Interestingly the rescaled Breit-Wigner approximation describes 
the POE data almost as well. Only for k = 1 we find minor deviations at intermediate times. This reaffirms, 
that the rescaled Breit-Wigner approximation is best suited for cases where the correlations are relatively weak. 
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Figure 9: As in Fig. ||, but for the cross correlation function divided by T^. Note, however, that in part (d), the 
C-function in the GOE case is exactly zero, and not plotted. 



Cross correlation functions 

Let us again think of a discrete stick spectrum where the levels are weighted by intensities (or partial widths). 
In order to calculate the C-function between two spectra related to different channels, it is reasonable to assume 
that the positions of the resonances are the same, but the corresponding intensities are uncorrelated. Then this 
cross C-function becomes: C{t) ^ 1 — 62(^)7 *-e. we obtain the full correlation hole |Q (see also the Eqs. (p3[)-(p7|) 
and the discussion following them). For the GOE case, this has been shown numerically in Sec. |]^. 

In Fig. |9| we show numerical and theoretical results for the cross C-function. In the GOE case we plot the 
numerical results together with the exact theory, Eq. (|3^), and we find again perfect agreement well within the 
statistical error (note that for k = 1, the cross correlation function vanishes identically). In the POE case we 
compare our numerical data with the rescaled Breit-Wigner approximation, Eq. (^7|). In the panels (a) and (b), 
which correspond to the weak coupling case, the theory agrees very well with the numerical data. For k = 0.3 (c), 
however, we find first systematic deviations. The rescaled Breit-Wigner approximation begins to overestimate 
the true C-function, which becomes even worse for k = I (d). 

Possibly the deviations for k > 0.3 may be due to the appearance of correlations when the resonances begin 
to overlap. The fact that the deviations are more pronounced here, than in the case of auto C-functions can 
be explained in a natural way by the higher sensitivity of the cross C-function to correlations. Moreover, it has 
been proven in Ref (though only for the one channel case) that the coupling to decay channels may indeed 
induce correlations even in the regime k < 1. 



In the GOE case, we know from Sec. 5.1 that the rescaled Breit-Wigner approximation reproduces the exact 
value of auto and cross C-functions at t — for arbitrary coupling strengths. This can also be proved for the 
POE in the one channel case |43 . Here, Fig. ^ and Fig. || give numerical evidence, that the same is true in 
the many channel case as well. Thus the relative size of the correlation hole seems to be independent of the 
coupling strengths, which turns the correlation hole into a persisting signature of chaos even in the limit of strong 
transmission Va : — > 1. This may be contrasted to the behavior in the energy domain, where the correlation 
hole for C[(Tfrt\ crfot ](w) gradually disappears with increasing number of channels or total transmission. Note, 
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Figure 10: The autocorrelation function for the elastic cross section divided by with the coupling strength 
K = 0.1. The numerical data are shown by circles (GOE) and crosses (POE). The theoretical curves are the 
rescaled Breit-Wigner approximation, Eq. (|3^), for the GOE (solid line) and the POE (dashed line). 



however, that in the time domain, the correlation hole becomes compressed to short times. One would then 
need very long spectra with many resonances to resolve the correlation hole. The generality of this invariance 
property and its implications will be the subject of future studies. 



7.2 Partial cross sections 

In this section we show results of Monte Carlo calculations for C-functions of partial cross sections. Though 
we obtained a closed formula even here (using the diagonal approximation for the partial cross sections and the 
rescaled Breit-Wigner approximation for the correlation function), we must be prepared for the approximation 
to deviate relatively early from the true result. The main advantage of partial cross sections, is their easy 
experimental accessibility. According to our approximation, the relative size of the correlation hole is given by: 

B (l + 6/Af)(l + 4/M) 

A (1 + 2/MY ^ ' 

(see Eq. ( |3l| ) and the list below). Note that for particular choices of the partial cross sections, this quotient may 
become even larger than one, which implies that the corresponding C-function for the GOE becomes negative 
at small times. 



Autocorrelation functions 

According to Eq. (|T]) we have to distinguish between the auto C-function of an clastic or inelastic cross section. 
In the elastic case, B/A — 3/35 so that the relative size Q of the correlation hole varies from Q = 3/35 « 0.0857 
for M ^ oo, and 9/35 w 0.257 for M = 2. In the inelastic case, B/A = 1/9 so that Q varies from Q = 1/9 « 0.111 
for M -> oo, and Q = 1/3 w 0.333 for M = 2. 

In Fig. |l^ (elastic case) and Fig. |l^ (inelastic case) we set the coupling parameter to k — 0.1. The numerical 
data is shown together with our approximate result, Eq. (pl|). In both figures, the approximation clearly deviates 
from the numerical result (which we may consider to be exact). The deviations are more pronounced in the 
GOE case, which is in line with earlier observations in the context of total cross sections. The deviations occur 
just in the region of the correlation hole, and thus diminish the difference between GOE and POE considerably. 
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t 

Figure 11: As Fig. but for inelastic cross sections. 



In Fig. |l^ it seems, that the theoretical values at t — are still quite accurate, so that the correlation hole could 
be detected there. However, the curves for the GOE and POE merge very rapidly (at t « 0.5). In Fig. |l^ we see 
that our approximation fails even at t = 0. The numerical curves for POE and GOE lie so close together, that 
it is practically impossible to identify the correlation hole. Similar investigations |8| gave essentially the same 
results. 

Cross correlation functions 

As the number of channels is M = 2, we need to consider only two different types of cross C-functions: the 
C- function of two elastic cross sections, and of an elastic and an inelastic cross section (sharing the same entrance 
channel). In the list below Eq. ( |3l] ) wc find that B/A = 1 in the first case, and B/A = 1/5 in the latter. Hence 
we may expect to obtain considerably larger correlation holes. For the case of two elastic cross sections the 
relative size Q of the correlation hole is Q = 1 for M — > cx) and Q = 3 for M = 2. This means, that the 
C-function in the GOE case becomes negative at small times. For the case of one elastic and one inelastic cross 
section, Q = 1/5 for M — ^ oo and Q = 3/5 for M = 2. This is the most natural constellation in experiments. 
Unfortunately the corresponding correlation hole is not very large. 



Figure 12 shows the cross C-function of two elastic cross sections divided by T^. In the four panels (a)-(d) 
the coupling strength is increased from k — 0.03 up to k = 1. The correlation hole is so big {Q — 3), that the 
GOE curve becomes indeed negative at small t, with the only exception of panel (d). The rescaled Breit-Wigner 
approximation works quite well in (a) and (b), begins to fail in (c), and fails completely in (d). Surprisingly 
the error becomes large in the POE case first, while the agreement in the GOE case is still reasonable. This 
contradicts to some extent, what we have found and partially explained in the case of total cross sections. Further 
studies are probably necessary to clarify this point. 

Figure |l^ shows the cross C-function of an elastic and an inelastic cross section. Not only is the correlation 
hole smaller, the deviations between our approximation and the numerical data are also larger than in the pre- 
vious case (Fig. |2|). In panel (b) for k = 0.1 we find already noticeable deviations in the POE case, and the 
correlation hole has practically disappeared for k = 0.3, in panel (c). In panel (d) finally, the numerical data 
practically coincide for the POE and the GOE case. It seems that the C-function in this case is negative for any 
value of t. Again this is an aspect, which deserves further studies. 

In this subsection, we saw that the validity of our approximation for partial cross sections is restricted to 
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Figure 12: The cross correlation function C[aii, (T22](i) divided by T^. The coupling strength is varied: k = 0.03 
(a), K — 0.1 (b), K = 0.3 (c), and k = 1 (d). The numerical data is shown with circles (GOE) and crosses (POE). 
The rescaled Breit-Wigner approximation ( ^ ) is plotted for the GOE (solid hne) and for the POE (dashed hne). 
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(b) ft = 0.1 



(d) K = 1 



Figure 13: As in Fig. [l|, but for the cross correlation function between an elastic and an inelastic cross section 
C[aii,ai2]{t) divided by T^. 
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coupling strengths that are typically much smaller than in the case of total cross sections. Moreover, we can 
no longer obtain the exact value of the C-functions a,t t — 0. However, we found here the largest effects of 
the correlation hole. In particular, the case of two elastic cross sections is very promising in view of a possible 
realization in an experiment. 



8 Summary 

By introducing an appropriate rescaling of the coupling strengths in the conventional theory of multichannel 
scattering for weak absorption, we have been able to extend its domain of validity far into the region of overlapping 
resonances. We have used this method (which we call the rescaled Breit-Wigner approximation) to calculate 
auto and cross correlation functions for total and partial cross sections both for the POE and the GOE case, i.e. 
in situations we consider as "integrable" and "chaotic." 

Starting from the VWZ integral which describes the GOE case, we have performed exact calculations 
for the correlation functions of total cross sections in the time domain. By comparison with these results and 
numerics, we determined the domain of validity of the rescaled Breit-Wigner approximation. It turned out, that 
it works particularly well in the case of weak correlations. 

We showed the differences between POE and GOE all the way up to the domain of strong transmission in 
all channels, and we found that autocorrelation functions are a poor means of determining the chaoticity or 
integrability of a scattering system in strong absorption, whereas certain cross correlation functions are quite 
sensitive. A particularly strong effect is seen in cross correlation functions of two different total cross sections. 
Unfortunately total cross sections are usually very hard to measure and, therefore, this observation may be of 
little practical value. As far as partial cross sections are concerned, we find that cross correlations of two elastic 
cross sections are also very sensitive, far more than either the ones between an elastic and an inelastic cross 
section or between two inelastic ones. 

For systems where two asymptotic channels are stable the cross correlation of two elastic cross sections is 
a very attractive means of analyzing the reaction. Also at least for intermediate coupling, i.e. overlapping 
resonances with moderate coupling in each channel we have a theoretical result with which to compare. This is 
particularly important for the integrable case where we know that nongeneric behavior might be quite common. 
First experiments on microwave billiards to check our results are on their way. 
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Average S-matrix 

The S-matrix given in Eq. (]^) may equivalently be expressed in terms of the so called "K- matrix" fl^ : 

As we will see below, the average S-matrix can be obtained easily from the average K-matrix. This offers a 
convenient way to define the openness of the system, i.e. the coupling strengths to the decay channels. 

In the eigenbasis of _ffi„t we get the following expression for the elements of the K-matrix (|6l|): 

1 ^ T/. T/. 

Ka,iE)^lY.^^ (62) 

i—l * 

where {si} are the eigenvalues of the closed system -ffi„t. The coupling matrix elements {Via} and the eigenvalues 
{si} are statistically independent in the ensembles considered (the GOE and the POE). Therefore, we get for 
the average of Kab[E): 

{Kab{E+)) = Sat y {V,a V,b) (^ j^_1^-q ) : (63) 
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which holds for any i. Here (. . .)p stands for the average over the level density p. Using that (V^a^fc) — 
5ij 5ab(y^a)-' which again holds for the scattering ensembles considered here (see Sec. 0), we get: 



(64) 



where we split the integral used to average over the level density, into its principle value part and the residue. 
Note, that the average K-matrix is diagonal. We define the center £^o = of the spectrum as that point where 
the principle value integral vanishes. Then we obtain: 

ITT 

(ifaa(0+)> - -- N{V^,) p{Q) . (65) 
Finally we define the "coupling parameters" K(i 8jS follows: 

'«a=i(ifaa(0+)) = ^ ^ N^) ' ^^^^ 

It can be shown that due to the analytic properties of the S-matrix its average is directly related to the 
average K-matrix: 

^^-iC^GBl)) 

^^(^^^ - l + i{K{E+)) ■ 

The average S-matrix is also diagonal. In the center of the spectrum the average K-matrix is purely imaginary, 
so that {S{E)) is real. Its elements are then given by: 

(SaaiO)) = . (68) 
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